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1. Introduction 

Lattice QCD simulations necessarily take place in a finite volume. This gives unphysical 
effects that must be compensated for. For example, using periodic boundary conditions removes the 
actual boundary in order to eliminate boundary effects but introduces conditions on the momenta 
available in the simulation. At a lattice size of 5 fm this means that the lowest lying non-zero 
momentum is 


p(l) «250MeV. (1.1) 

Depending on the process under study this can be considered large, and furthermore there are 
situations where being able to continuously vary momenta is advantageous. Examples include 
simulations for the charge radius of the pion |[T]], hadronic vacuum polarization for /r g-2 0 and Kj-i 
decays close to zero momentum which is relevant for determination of Vus 1^- Twisted boundary 
conditions, which are introduced in section ^ is a way of continuously varying momenta in lattice 
QCD [| I i. 

In this talk we investigate some effects of twisted boundary conditions using continuum SU(3) 
Chiral Perturbation Theory in finite volume at one loop ^ |^. In particular we show how 
Ward identities are satisfied in a nontrivial way and how this relates to meson masses. Full results 
for all calculations will not be presented since they are already available in [[TT[]. We focus instead 
on the principal differences between the case of periodic boundary conditions and twisted boundary 
conditions using results for 71+ to illustrate. See also [ |I^ for discussions on some of the issues 
covered here. 


2. Twisted boundary conditions 


Twisted boundary conditions are defined through the relation 

+L) = 0(x')exp(/0'). (2.1) 


where we have assumed that the field is confined to a cube where all sides are of equal length L 
and with i running over spatial dimensions. Expanding both sides in Fourier modes gives 


f' n 97r 

dp'^{p')e.x^{ip\x' +L)) = j dp'(j){p')exp{i{p‘x‘ + 6)) p' = p'{n') = -j-n' -f—,ne (Z). 

( 2 . 2 ) 


By varying the twist angle, 6, it is possible to vary the momentum of the field continuously. There 
are two more important observations we want to point out. First, twisted boundary conditions break 
the cubic symmetry of the lattice which means that e.g. form-factors are not as constrained. Second, 
available momentum states are no longer symmetric around zero. This has the consequence that 
tadpole integrals are no longer zero and, more generally, that there are new terms in loop integrals 
which reflects the loss of cubic symmetry. These two related observations were the starting point 
of our calculations. 
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2.1 Twist for quarks 

The Dirac term in the QCD lagrangian is 


= q{ip — m)q (2.3) 

where q and q are Dirac spinors and also vectors in flavour space, D is the covariant derivative in 
QCD and m is the flavour diagonal mass matrix. We now define twisted boundary conditions for 
each quark flavour separately as 

qf{x' +L)= exp(ie])^/(x'). (2.4) 

This definition makes it clear that each flavour and direction can have its own twist angle To 
make contact with methods developed for periodic boundary conditions it is practical to undo the 
twist of the quarks using a field redefinition [^]. This redefinition can be seen as a local U{n)v 
flavour transformation for n flavours, using only the Cartan subalgebra. It is however not to be 
confused with a gauge transformation in color space. Collecting the twist angles for the flavour / 
in a three vector 6f we define 

qf{x) = qf{x)exp {-Wf ■ ^ . (2.5) 

The field q satisfies periodic boundary conditions 

q{x' +L) = q{x‘). (2.6) 

The lagrangian for the field q is then 

.if = q{ip + — ni)q (2.7) 

where we’ve collected the twist angles in a flavour diagonal matrix 0^ where the entries are four 
vectors 6^ = (0,6/). We can now interpret the term qif)/Lq as coming from a constant background 
field. Note that it interacts differently with the different flavours since it is proportional to the twist 
angles. Note that the quarks, q, have momenta k = ^n. 

2.2 Twist for mesons 

A meson with the flavour structure q'q has the twist angle 

^qiq{x'+L) = exp(/(e' - 0',))0^/^(x') = exp(/e|,^)(/)^/^(x'). (2.8) 

We will use the three- and four-vector notation for meson twists as we did above for the quarks 
and we will label the mesons with .... Note that even though the flavour diagonal mesons 
have twist angle zero, they are affected through mixing in the cases where twisted boundary con¬ 
ditions break isospin symmetry. The broken isospin symmetry affects the decay constants at NLO 
while the effect for the masses is NNLO [p^]. For charged particles charge conjugation is broken. 
This happens since particle and anti-particle have opposite twist angles. A particle with spatial 
momentum p corresponds to an anti-particle with spatial momentum —p. 
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2.3 Calculating Feynman diagrams with twisted boundary conditions 

We will now introduce the integral notation used in the rest of the talk, for a more complete 
treatment see Bin HI- We label tadpole integrals with A according to 




d‘^k 

{IkY {k^—iffiY 


(2.9) 


and two propagator integrals with B according to 


B 






mi ,m 


1 5 "^2 5 



d‘^k 

(IzY {k^ — mjY' (id — pY 



( 2 . 10 ) 


When spatial momenta are quantized, the spatial part of loop integrals are turned into sums 
over spatial momenta as 

fr fy 
J {lltY y (2l)''-3 ■ ' ' 

' ^ ^ ^ ^ ^ n(^I?,k=(2nn+e^)lL 

In it is explained how all propagators in a Feynman diagram end up with allowed momenta. 

The sum of course modifies the numerics of each integral but also explicitly breaks cubic 
invariance. This means that the finite volume parts of loop integrals, labeled by and , do not 
decompose into Lorentz structures in the usual way. Most notably 


A^{m\,n) 


d^k 

{InY {k^ — tn^Y 


( 2 . 12 ) 


so that the odd tadpole integral is non-zero even though there is no Lorentz structure available. This 
happens due to the fact that momenta are not symmetric around zero so that the usual cancellation 
does not occur. Other decompositions that we use in our results are, using 1,1) = 

5{.M.MV}(^^^2), 

{rn\,rnl) = q^B\{rn\,rnl,q) +bI^ {rn\,rnl,q) 

B^'^^ {m\,ml) = q^ q'’BlYm\,mYq) + {ml,rnl,q) + bI^''{ m\,ml,q) 

(2.13) 


For the untwisted case the integral b'^^ = 0 while B^^'' is zero in infinite volume. Note that the 
integrals on the right hand side depend on the components of q, not only on q^ as would be the case 
in infinite volume. Note that all integrals containing twisted mesons depend on twist angles in the 
finite volume results. It is therefore important to specify whether there is a 71+ or a 7i^ in a loop 
and similarly for other particle anti-particle pairs. 

When presenting our results we will use the notation 

A^X=X{V)-X{oo) (2.14) 

where X is the object under discussion. We call A^X the finite volume correction. 
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3. Masses 

In order to discuss masses, we first need to define whaf we mean by mass. Our definifion is 

m^=E^-f (3.1) 

for fixed spatial momentum p. This will result in a momentum dependent mass m = m{p) due to 
the broken cubic symmehy. Our definition should be conhasted with the definition [ |l2| ] 

m^ = E'^-{p + Kf, (3.2) 

where is a NLO correction to p which removes the momentum dependence. The two definitions 
differ simply by moving around parts of the loop corrections. Note that the shift in momentum is 
not universal [|TT|]. With this we mean that, if we define the decay constant as 

(OIA" \M{p)) = iVlEM (p^ +K’^) , (3.3) 

in order to absorb the momentum dependence of the decay constant we would need to renormalize 
the momentum differently, that is K' / K. Note also that it is not possible to find any K such thaf 
the elechomagnetic form-factor can be written as a single /+ with rescaled momentum. 

Now that we have defined mass, we quofe the one loop result for the pion mass to illushate the 
differences that come up when using twisted boundary conditions; 



The factor in the first row inhoduces momentum dependence to the mass. All integrals in the 
first row vanish as the twist angles go to zero. Note that the mass shifts for a 71+ with momenta p 
and a with momenta —p are the same which is consistent with that charge conjugation should 
be defined wifh a minus sign for momenfa. The correcfions for flavour neuhal mesons remain 
momentum independent. 

4. Electromagentic form-factor and Ward identities 

We define the electromagnetic form-factor by 

{M' {p)\fllW{p)) =flMM'^l 

= flMM' + {Pll + P'p) + flMM'-{Pll — Pn)+hMM'^ ■ (4-1) 

where M and M' are mesons and is the electromagnetic current. We split up the form-factor into 
three parts with different momentum dependence. The split is not unique but comes naturally when 
calculating the finite volume corrections. As long as we only look at the elechomagnetic current. 
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the Ward identity ensures that fiMM'- and hjMM'^i go to zero when the twist angles go to zero. For 
the charged pions the one loop finite volume corrections are 

^^fem 7 t±+ = ^ +H'^{m\+,m\-,q)) 

v\ 

^^femn+{n-)- = -( 2 ^^ > " 4 -> <?) + ^Iv (^1+, ml-, q)) 

^0 


,Vi _ _j_ 

^ ^emn^Li 

^0 


2A)i (m2 +) + a|^ (m|+) - (m|o) 


*? ^2)1 i^7C+ 5 ^71^ ) ^) “h 2 ^2^ (^A'+ > 

TiP + p'V (2^23^V {^1 +, ml- , ^) + B23mV ('”i+ ’ > <?)) 


(4.2) 


where, apart from the integrals defined in secfion [2.3| , we have used a generalizafion of found 

in|H,|il 




(4.3) 


The Ward identify 

{Pii — p'ii)flMM'n = 0 (4.4) 

is still satisfied. This does fake some exfra care as we show below. Expanding fjMM'^i we gel 

0 = qn {M'{p)\jl\M{p)) = fiMM'+{P^ -P'^) + fiMM'-q^ + hiMM'^qii (4.5) 

where q^ = (p^ — p'p)- In the infinite volume case for M = M' we always have (p^ — p'^) = 0 on 
shell fo all orders. In fhe case of Iwisfed boundary conditions fhis is no longer frue. Inslead we 
have, using 71+ as an example. 


2 /2 
P -p 


NLO 


ph — p'h 
P’2 


[-2A)((m2 + ) -Al{ml+) +A);(m^)] 




\V t 2 


(4.6) 


which gives a non-zero confribufion when multiplied wifh fhe leading order pari of fiMM'+- This 
confribufion cancels parls of hjMM'^q^i- Thereminaing terms cancel befween//M;vf /_(72 and hiMM’yLqii 
using some infegral idenfilies [pT|]. 

An issue which is nol solved by using twisled boundary condilions is lhal for M = M' Ihere is 
no way lo continuously vary fhe momenlum in fhe currenl since p‘ — p‘f = — n‘r). In order 

lo address Ibis case if has been proposed in [1^] lhal broken isospin symmelry should be used lo 
relate malrix elements through 


{n+{p')\uY^u\n+{p)) = {n°{p')\dY^u\n+{p)). 

Note that in this way contributions from the sy^s part of the current are missed. 
We now split the form-factor as 

fn = (!+/+ +AV+) {p + p )ti+^^f-qti +A'^/7^ = 


(4.7) 


(4.8) 
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Figure 1: Left: = 0 Right: = 1. One loop corrections to /+. Plotted are the infinite volume correction, 
ff{q^), the finite volume correction, A'^/+, and the full finite volume correction, A^f^ = A'^/+(p + p')^ + 
A^ -\-A^h^. The one loop finite volume corrections are of similar size as the infinite volume one loop 
corrections. 

where we have separated out the lowest order value of 1. In Fig. [I] we plot the time and space 
components of the Afferent parts of /+ as well as the full finite volume correction in a volume with 
size = 2. In order to look only at the loop contribution we have set Lg = 0. The renormalization 
scale is set to jj. = 0.77 GeV and we have used the following input constants 

m;r = 139.5 MeV, m/f = 495 MeV, = Fj, = 92.2MeY, mj,L = 2. (4.9) 

The momentum exchange is varied by changing the twist angle on the u-quark. From Fig. it is 
clear that the finite volume loop correction is of a similar size as the infinite volume loop correction. 


5. Decay constants 


Similar to the case of the electromagnetic form-factor, the decay constant also has a different 
momentum dependence with twisted boundary conditions than with periodic boundary conditions. 
We define fhe (axial-vecfor) decay consfanf in finife volume as 


(OIa" \M{p)) = iV2FMPy, + iV2F^^ , (5.1) 

where M{p) is a meson and wifh fhe corresponding axial-vecfor decay consfanf = qy^y$ (A^/ V2)q. 
We also define a pseudo-scalar decay consfanf 


(0|P“|M(p)) 


Gm 

V2 


(5.2) 


where = qiy^ {2.^/y/2)q is fhe pseudo-scalar current corresponding to the meson M. For flavour 
charged mesons composed of qq' the two matrix elements satisfy the Ward identity 


{Q\A^\M{p)) = {m, + mq,){Q\P^\M{p)) (5.3) 
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Figure 2: The relative size of the finite volume correction to the decay constant of plotted against m-^L 
for various twist angles on the u-quark. One the left we have plotted if'F m/Fji and on the right we have 
plotted F}^^/{F„m„). 


which leads to 

P^Fm + = -{mq + mq')GM ■ (5.4) 

When checking this Ward identity, it is important to use the correct momentum dependent mass. 
The explicit expressions for decay are given by 


rfip = ±4 (2A);(mJ,) +Aj;{m|*) , 

A''G;;* = ^(iA'-(4,) + iA''(4.) + iA''(;„i,) + lA''(™J)) . (5.5) 


The finite volume corrections to 7r+ are plotted in Fig. ^ 

As a final note on decay consfanfs we note fhaf due fo fhe isospin breaking effecf of fwisfed 
boundary condifions fhere is mixing in fhe decay of fhe flavour neufral mesons. We will nol go info 
fhe defails here buf refer fo [|Tl|] for fhe one loop resulfs and | lA] for how fo freaf mixing. 


6. Conclusions and outlook 

In Ibis falk we have presenfed parts of our work on fwisfed boundary condifions using SU (3) 
ChPT. The sfarfing poinf is fhe observafion fhaf fhe odd fadpole integral is no longer zero when 
fwisfed boundary condifions are used. We note also fhaf fhe broken lallice symmefry means fhaf 
exfra care is needed when defining quanfifies such as form-factors. 

More specifically, we have given an infroducfion fo fwisfed boundary condifions and how 
compufafions are affected. We have also presenfed resulfs for mass, elecfromagnefic torm-facfor 
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and decay constant for Finally we illustrated how the relevant Ward identities are satisfied. The 
procedure for checking Ward identities can be summarized as be careful and do things correctly! 

The size of the effects of twisted boundary conditions can be considered small in many cases, 
but depending on lattice size and the precision desired they can be non-negligible. 

We are now working on the one loop corrections to Klj, decay in partially twisted partially 
quenched staggered chiral perturbation theory, relevant for determining the Vus mixing angle [[I^. 

Acknowledgments 

We thank the organizers for a very pleasant and well-run conference. This work is supported 
in part by the Swedish Research Council grants 621-2011-5080 and 621-2013-4287. 

References 

[1] B. B. Brandt, et al.. The pion vector fonn factor from lattice QCD and NNLO chiral perturbation 
theory, JHEP 1311 (2013) 034 [arXiv: 1306.2916 [hep-lat]]. 

[2] C. Aubin, et al, “Hadronic vacuum polarization with twisted boundary conditions,” Phys. Rev. D 88 
(2013) 7, 074505 [arXiv;1307.4701 [hep-lat]]. 

[3] P. A. Boyle et al, “The kaon semileptonic form factor in Nf=2H-l domain wall lattice QCD with 
physical light quark masses,” [arXiv: 1504.01692 [hep-lat]]. 

[4] P. F. Bedaque, Aharonov-Bohm effect and nucleon-nucleon phase shifts on the lattice, Phys. Lett. B 
593 (2004) 82 [arXiv:nucl-th/0402051] 

[5] C. T. Sachrajda and G. Villadoro, Twisted boundary conditions in lattice simulations, Phys. Lett. B 
609 (2005) 73 [arXiv:hep-lat/0411033]. 

[6] G.M. de Divitiis, R. Petronzio, N. Tantalo, On the discretization of physical momenta in lattice QCD, 
Phys. Lett. B, 595 (2004) 408 [arXiv:hep-lat/0405002]. 

[7] S. Weinberg, Phenomenological Lagrangians, Physica A 96 (1979) 327. 

[8] J. Gasser and H. Leutwyler, Chiral Perturbation Theory To One Loop, Annals Phys. 158 (1984) 142; 

[9] J. Gasser and H. Leutwyler, Chiral Perturbation Theory: Expansions In The Mass Of The Strange 
Quark, Nucl. Phys. B 250 (1985) 465. 

[10] J. Gasser and H. Leutwyler, Spontaneously Broken Symmetries: Effective Lagrangians At Finite 
Volume, Nucl. Phys. B 307 (1988) 763. 

[11] J. Bijnens and J. Relefors, Masses, Decay Constants and Electromagnetic Form-factors with Twisted 
Boundary Conditions, JHEP 1405 (2014) 015 [arXiv: 1402.1385 [hep-lat]]. 

[12] F. -J. Jiang and B. C. Tiburzi, Phys. Lett. B 645 (2007) 314 [hep-lat/0610103]. 

[13] G. Amoros, J. Bijnens and P. Talavera, QCD isospin breaking in meson masses, decay constants and 
quark mass ratios, Nucl. Phys. B 602 (2001) 87 [arXiv:hep-ph/0101127]. 

[14] D. Becirevic and G. Villadoro, Impact of the finite volume effects on the chiral behavior of fK and Bk, 
Phys. Rev. D 69 (2004) 054010 [arXiv:hep-lat/0311028]. 

[15] J. Bijnens, E. Bostrom and T. A. Lahde, Two-loop Sunset Integrals at Finite Volume, JHEP 1401 
(2014)019 [arXiv: 1311.3531 [hep-lat]]. 


9 



Masses, decay constants and em form-factors with twisted boundary conditions 


Johan Relefors 


[16] J. Gasser and H. Leutwyler, Low-Energy Expansion of Meson Form-Factors, Nucl. Phys. B 250 
(1985)517. 

[17] J. Bijnens and R Talavera, Pion and kaon electromagnetic form factors, JHEP 0203 (2002) 046 
[arXiv:hep-ph/0203049]. 

[18] J. M. Flynn et al. [UKQCD Collaboration], A Numerical study of partially twisted boundary 
conditions, Phys. Lett. B 632 (2006) 313 [arXiv;hep-lat/0506016]. 

[19] C. Bernard, J. Bijnens, E. Gamiz and J. Relefors, work in progress. 


10 




